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SUMMARY 


A theoretical method is presented for the computation of the flow field 
about an axlsymmetric body operating in a viscous incompressible fluid. This 
approach combines a smoothing routine, a potential flow method based on a 
surface source distribution, and a f Inlt e-dlf f er ence boundary- layer method to 
accomplish the analysis. An empirical method used for modeling separated 
flow is shown to work reasonably well for cases of extreme flow separation. 
Results obtained by this method are presented which show very good agreement 
with experimental data. Suggestions are made for extending this method both 
to Include a better model for separated flow and to calculate the "viscous" 
flow about axlsymmetric bodies at angle of attack. A detailed instruction 
manual for input Ing data to the computer program is given in Appendix A. 
Appendix B contains the necessary information to place this program on to a 
computer. This appendix also contains a complete description of output para- 
meters from the computer program, as well as basic flow charts of some of the 
major subroutines. Appendix C contains a complete listing of the computer 
program for operation on either a CDC or an IBM computer. 
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A TIffiORETICAL METHOD 
FOR THE ANALYSIS AND DESIGN 
OF AXISYMMETRIC BODIES 


by T. D. Beatty 

McDonnell Douglas Aircraft Corporation 

One of the ultimate goals in aerodynamics is the achievement of the 
ability to obtain the real fluid flow field about an arbitrary three- 
dimensional configuration by theoretical calculation rather than by resorting 
to expensive and time consuming wind tunnel tests. The exact treatment of 
this problem requires the solution of the full Navier-Stokes equations, which 
is currently not practical. However, a good approximation to this real flow 
can be obtained by displacing the surface boundaries of the original body 
to account for viscosity as shown by Thwaltes in Reference (1), 

This technique of displacing the boundary surface to obtain a viscous 
solution has been used in two-dimensional flows quite successfully, as shown 
in References 2, 3, and 4, The extension of this approach to three-dimensional 
flow requires that appropriate computational routines be available to calculate 
the potential and viscous flow parameters. A potential flow routine which can 
calculate the flow about arbitrary three-dimensional bodies is available 
(Reference 5), although the con 5 >arable three-dimensional boimdary layer method 
is not currently in the state of the art. At the present time, the general 
three-dimensional problem cannot be solved. However, both an axiS 3 anmetrlc 
potential flow method and an axisymmetric boundary layer method which can 
calculate the Invlscid and viscous flow field about a body of revolution at 
zero degrees angle of attack are currently available. 

Because of its simple nature and its common appearance in fluid dynamics, 
it was decided that a body of revolution would be a good starting point for 
the development of a three-dimensional method for calculating invlscid and 
viscous flow fields. 

The axisymmetric potential flow routine (References 6 & 7), used in the 
present method was developed at the Douglas Aircraft Conq>any under the guidance 
of A. M. 0, Smith and has proven over the years to be an extremely versatile 



and accurate method, as well as the only purely axisymmetric potential flow 
method, generally available in industry today. This method has been well 
disseminated throughout industry; only a brief discussion will, therefore, 
be presented in a following section. 

The boundary layer method presented in this report (Reference 8) is a 
finite difference technique which uses an eddy- viscosity concept to replace 
the Reynolds shear stress term. Since this method is relatively new and has 
been modified extensively since Reference 8 was reported, a detailed 
description will be presented. 

The capability of the present method to determine the viscous flow 
about axisymmetric bodies is shown by correlations between the calculated 
results and experimental data. 

Recommendations are presented for extending the present method to the 
calculation of the flow field about axisymmetric bodies at angle of attack. 
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DEFINITION OF SIMBOLS 


A 

A+ 


D 

f 

G 

H 

k 

L 

Z 

P+ 


R 

X 


r 


T 


t 


U 


00 


U 

T 

U 


Damping length or frontal area, wherever applicable 

Damping constants 

Total skin friction coefficient 

Pressure coefficient 

Chord 

Local skin friction coefficient 
Maximum diameter 
Dimensionless stream function 
/Spot formation parameter 
Shape factor, 0/6 
Mixing-length constant 

Power to determine 2-D or axisymmetric flow 
Reference body length 
Mixing length 

Pressure gradient parameter 
Chord Reynolds number, u^c/v 
Diameter Reynolds number, 

Local Reynolds number. Ugx/v 
Momentum thickness Reynolds number, V^B/v 
Radial distance from axis of revolution 
Local radius of body of revolution 
Absolute temperature, ®K or "ET. 

Transverse curvature term 
Free stream velocity 
Friction velocity, 

X component of velocity 
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Velocity at edge of boundary layer 

V y-component of velocity 

X Distance along surface measured from leading edge or from 

stagnation point 

y Distance normal to the surface of the body 

a Angle between normal to the surface y and the radius r 

“ Constant in outer eddy viscosity equation 

6 Dimensionless velocity - gradient term, B = (2g/Ug) (du^/dg) 

Transitional parameter 
6 Boundary layer thickness 

6 Boundary layer displacement thickness 

e Eddy viscosity 

4 * 

e Ratio of eddy viscosity to kinematic viscosity, e/v 

n Transformed y~coordinate 

0 Momentum Thickness 

M Dynamic viscosity 

V Kinematic viscosity 

C Transformed x~coordinate 

p Density 

T Shear stress 

SUBSCRIPTS 

c Switching point between the inner and outer eddy viscosity formulas 

e Outer edge of boundary layer 

1 Inner region 

A Laminar 
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o Outer Region 

t Turbulent 

Tr Transition 

w Wall 

■=- Free-stream conditions 

Primes denote differentiation with respect to n. 
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TECmiCAL DISCUSSION 

Geometry Definition 

The geometry input to the Douglas Neumann Potential Flow Program must 
satisfy two primary requirements! the coordinates must be distributed 
properly and the surface curvature must be smooth. These requirements are 
easily achieved on an analytical body shape, since the input coordinates may 
be calculated exactly for any prescribed distribution. However, some method 
of determining accurate input coordinates for an arbitrary axisymmetric body 
is necessary, since the body may not always be amenable to exact analytical 
definition. The approach adopted in the following method is to assume that 
the coordinates are input in the proper distribution about the body, but that 
they are not necessarily smooth. These two requirements will be discussed 
in some detail in the following sections. 

Point distribution . - In order to obtain a high degree of accuracy in 
defining a pressure distribution when using the Douglas Neumann Potential 
Flow program, surface coordinates should be concentrated in regions of high 
surface curvature where rapid changes in the surface pressures would be 
expected. Since the total number of points per body is fixed, the distri- 
bution of these points about the body contour becomes extremely important. 

The Neumann program uses the input coordinates to create linear segments 
between points , thus approximating the body by a series of Frustums of Cones. 
The basis distribution required is then quite simple! more points and thus 
smaller segment sizes in regions of high curvature and less points and thus 
larger segment sizes in the other areas of the body. The basic guidelines 
to follow to insure proper point distribution are simply that the surface 
lengths of adjacent elements should not change by more than twenty to thirty 
percent and the maximum length of any segment should not exceed either five 
percent of the body chord or fifty percent of the local body thickness. 

Smoothness of input coordinates . - The Douglas Neumann program, or any similar 
potential flow method, is sensitive to the derivative of the surface slopes, or 
the curvature of the surface. The surface defined by the input coordinates must 
therefore have smooth first and second derivatives. The approach used in the 
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present method to smooth these coordinates, is a five point smoothing routine, 
which assumes that the input coordinates are smooth and continuous to graphical 
accuracy, l.e. , points are chosen from a small graph (approximately a 10 inch 
chord) . The output points from this routine will be moved very slightly to 
smooth the derivatives, but this movement will be negligible as far as the 
body shape is concerned. The equations used to accomplish this smoothing are 
as follows : 


- A I 

“ 16 j" *J-2 

where x^ and 

coordinates 

and It, and 'y , 
3 3 


+ + lOXj + 4xj^^ - 


+ + lOy^ + - y^^^j 

y^ are the unsmoothed input 


are the smoothed coordinates, 


(la) 

(lb) 


Potential Flow Method 


ihe Douglas Neumann method, (References 6 and 7) is very general in 
that it can calculate the potential flow about virtually any body. There is 
no restriction, for example, to slender bodies; in fact, the ''body” In 
question need not be a single body but may be an ensemble of bodies. In 
principle, the calculated solution may be made as accurate as desired by 
suitably refining the numerical procedure; accordingly, the so-called Neumann 
method is designated an exact method in this sense. 

The Neumann method is based on the use of a distribution of source 
density over the body surface. Applying the condition of zero normal velocity 
on the body surface yields an Integral equation for^ the source distribution. 
Specifically, the equation is a Fredholm Integral equation of the second kind 
over the body surface. Once this has been solved for the source distribution, 
all flow quantities of interest, i.e., velocity, pressure, etc. , can be calcu- 
lated fay rapid straightforward procedures. To implement this method on a 
computer, the body surface is approximated by a large number of small surface 
segments, over each of which the source density Is assumed constant. The 
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integral equation is replaced by a set of linear algebraic equations for the 
values of the source density on the segments. Input to the computer program 
consists of the coordinates of a set of points defining the body surface; 
these points are then used to determine the surface segments for approximating 
the body. There is no assumption made that the body can be analytically 
represented. 

The usefulness of potential flow with its neglect of viscosity and 
compressibility is due to the fact that it is a good approximation to real 
flow under a wide variety of circumstances. With regard to viscosity, the 
program obtains useful results except in regions of catastrophic separation. 

To verify the usefulness of potential flow as a predictor of real flow, results 
calculated by the Neumann program have been compared with experimental data. 
Several collections of comparisons have been made. Reference 9 was a very 
complete collection but is now rather old. Reference 10 is a more recent 
collection that shows a smaller number of comparisons. In the calculation 
of the viscous flow about axisymmetric bodies it is necessary to add the 
boundary layer displacement thickness to the body as will be shown in a sub- 
sequent section. This results in an "open" trailing edge body. This "open" 
body can be evaluated by the Neumann program without any difficulty even 
though the boundary surface does not close. Reference 11 presents an 
explanation of this phenomenon which proceeds as follows : for a closed 

body the integral of the source density over the body is zero; for an "open 
trailing edge body, this integral is not zero, and a streamtube leaves the 
trailing edge of the open body which proceeds do\-mstream and approaches 
infinity parallel to u^ as a constant cross section streamtube. Thus, the 
flow that is calculated may be thought of as that about a semi- infinite body 
consisting of the open body and an extension defined by this streamtube. The 
shape of the extension is unknown but is presumably unique, having both zero 
normal velocity and zero source density. 

The potential flow program has many useful options available which do 
not pertain directly to the present development. The details of these options 
are described in References 12 through 17. 
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Boundary Layer Method 


Basic boundary layer equations . - The calculation of the viscous flow over an 
axlsyxranetrlc body Involves the solution of the laminar and turbulent flow 
equations. For laminar flows, the problem Is strictly mathematical because 
the governing differential equations can be written exactly. For turbulent 
flows on the other hand, an exact solution of the governing equations Is not 
possible. Consequently, In order to proceed at all, one must rely on a cer- 
tain degree of empiricism. In the past, most of the work In this area has 
concentrated on so-called momentum and/or energy Integral methods as a means 
of evaluating the viscous flow parameters. Thus, the exact mathematical 
solution to the problems of the turbulent flow was bypassed, leading to fast 
and simple methods with varying degrees of accuracy. These methods usually 
rely quite heavily on empirical correlations and generally are restricted to 
a limited range of flow conditions. 

The Douglas Boundary-Layer Method (Reference 8), eliminates many of 
the disadvantages of the Integral methods by proceeding to solve the full 
partial-differential equations governing the flow, thereby, being classified 
as a differential method . For two-dimensional and axlsymmetric Incompres- 
sible flows, turbulent boundary-layer equations contain terms involving time 
means of fluctuating velocity components known as Reynolds stress terms. At 
present the exact relationship between these terms and the mean velocity 
distribution in the boundary layer still remains unknown. In the present 
method, a relation based on the eddy-viscosity concept is used giving highly 
satisfactory results for a variety of flow conditions. 

If the normal-stress terms are neglected, the Incompressible turbulent 
boundary-layer equations for two-dimensional and axlsymmetric flows can be 
written as In Reference 8; 

Continuity 

^ P"] 'h P’] ■ ° 
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Momentum 


where 


u 


9u 

8x 


+ V 


3u 

3y 


u 

e 



»L Jl. _L 

Pco 9y 



with 





(For laminar flow only) 


( 3 ) 


(4) 


T “ - p u’v' (Additional term due to turbulent flow) 

t 

and 

u' v' “ Reynolds shear stress term 

k ■ 0 for two-dimensional flow 
k “ 1 for axlsymmetric flow 

The basic notation and coordinate scheme are shown in Figure 1, where U 
is a reference velocity and Ug(x) is the velocity Just outside the boundary 
layer. The coordinate system is a curvilinear one in which x is the 
distance along the surface measured from the stagnation point or leading edge, 
and y is measured normal to the surface. Within the boundary layer, the 
velocity components in the x- and y-directions are u and v, respectively. 
The body radius is r^. 


The boundary conditions for equation (3) are 

u(x,0) ” 0 
v(x,0) - 0 


u (x) 
e 


(5a) 

(5b) 

(5c) 


Llm u(x,y) 

y-xo 

Before equations (2) and (3) can be solved, they must be transformed to a 
coordinate system which removes the singularity at x = 0 and stretches the 
coordinate normal to the flow direction. First, these equations are placed 
in an almost two-dimensional fonn by the Probsteln-Elliott transformation 
(Reference 18) : 


dx 


dy 




dx 


dy 


( 6 ) 


(7) 
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where (x) is the body radius and r(x,y) is a radius which accounts 
for the transverse curvature effect which will be subsequently discussed . 

A stream function Y is defined that satisfies the continuity equation (2); 
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9y 


k 


r u 


H 

9x 


k 

r V 



( 8 ) 


The resulting equations are transformed by the Levy-Lees transformation 
(Reference 19) in order to remove the singularity at x = 0 and stretch the 
coordinates in the x and y directions. The Levy-Lees transformations are: 


00 00 0 
Pco'^co 

dn “ ' i dy 

(2KP 


(9a) 

(9b) 


A dimensionless stream function, f, is Introduced which is related to T 
as follows: 


T - (20* f(5,n) 


( 10 ) 


Combining the Levy-Lees and the Probsteln-Elliott transformations given above 
we have , ._2k 



J 

dx 

(11a) 

Poo'^e 

dn » —V 
(20* 

rr(x.y)p 

dy 

(11b) 


Introducing an eddy viscosity term to account for the Reynolds shear stress 
terms. 


u'v 

E , 


9u 

3y 


e = 


U 


( 12 ) 


and a transverse curvature term t along with a pressure parameter term 

O _ 2i ^ 


Ue dC 


(13) 


The momentum equation (3) then becomes, with f* = u/Ug, 


(l+e+) f"j + ff" + e^l-(f’) ^ = 2cj^f’ II' - f" llj (14) 
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The boundary conditions given by equation (5) become 


f(5,0) « f„ « 0 

(15a) 

f'(C.O) - 0 

(15b) 

Llm f ' (5 , n) = 1 

(15c) 




The momentum equation is then solved by a very efficient numerical scheme 
developed by Keller, (Reference 20) and applied to boundary layer calculations 
by Cebeci and Keller, (References 21 and 22). 


Eddy viscosity equations . - The eddy viscosity concept is used to relate the 
time-mean fluctuating velocities to a mean velocity distribution as given 
in equation (12) 



( 12 ) 


A two-layer model of the eddy viscosity within the boundary layer will be 
used as shown in figure 2. 


In the inner region of the boundary layer an eddy viscosity model, 
based on Prandtl's mixing-length theory, is used: 

« £2 |~ I (16) 

where £, the mixing length is given by 

£ - Kj^ Y (17) 


A modified expression for £ has been developed by Van Driest (Reference 23) 
to account for the viscous sublayer close to the wall. This modification is 

-(Y/A)-| 


where A is given by 


V [i- 


J 


(18) 


(19) 
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and 


A+ - 26.0 


= jl - 11.8 P+ 


V u 

00 e 


S' 


du 
e 

dC 





Now for axisynnnetric flows the 


X, = .4ro 


value of £ is replaced by 



( 20 ) 

( 21 ) 

( 22 ) 


(23) 


(24) 


which is developed in reference 24. If transverse curvature effects are 
desired then 


r 

r 

o 


r + Y cos a 
o 



1 


+ 


Y 

— cos a 
^o 


= 1 + t 


(25) 


where t 


Y 

— cos a 
r 

o 


then £ becomes 


£ 


.4r In (1+t) 
o 


- ^ ln(l+t), 

1 - e 


(26) 


The eddy viscosity in the outer region of the boundary layer is given by 


e 

o 


a u. 


(27) 


where 


is the boundary layer displacement thickness defined by 
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dn 


( 28 ) 



which in the transformed plane becomes 


n 


L 


( 25 ) 
P, 


^ / 


(1-f ') (1+t) dh 


(29) 


where 


1+t 



2L cos a 




(30) 


This relationship for Eq Is the same for two-dimensional or axisymmetric 
flows as shown in Reference 24 < 


Low Reynolds number effects . - The calculation of turbulent boundary layers 
about two-dimensional and axisymmetric bodies must often be done at low 
Reynolds number, i.e. , momentum thickness Reynolds number, Rg, less than 
6000. Most of the boundary layer methods including the one presented above 
are based on empirical data which were obtained at high Reynolds numbers. A 
correction term to account for low Reynolds numbers which was developed by 
Cebecl (Reference 25) based on prior work by Coles (Reference 26) is, there- 
fore, applied to the outer eddy viscosity by varying the a in equation (27) 
with R- in the following manner. 

O 

if R^ < 425 then a - (.0168) (1.55) (31a) 

if Rg > 6000 then a ■ .0168 (31b) 

if 425 < Rg < 6000 then a - .0168 [j“] where 


n 


.55 


-. 243/7 

1 - e 


.298 y1 


and Y 



(31c) 


Transverse curvature . - In developing the axisymmetric boundary layer 
equations a radius term is Introduced as shown in equations (2) and (3). If 
the assumption is made that the body radius is very large compared to the 
boundary layer thickness then the radii in equations (2) and (3) reduce to 
the local body radius r^ and the effect of the transverse (i.e. , 
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circumferential) curvature In the momentum equation is neglected. 
however, the body radius is small compared to the boundary layer thickness 

then the effect of the transverse curvature cannot be ignored and t must 
be a function of the distance into the boundary layer, y. The relationship 
between y, r^, and r is given by; 

r “ r + y cos a (32) 

As observed in figure 3, a Is simply the surface slope In the longitudinal 
direction, l.e. , 

tan « (33) 

For slender cylinders where a ■ 0“, 

r " r + Y (34) 

o 

The inclusion of the transverse curvature terms in the boundary layer 
equations is shown in References 24 and 27 to substantially Improve the 
accuracy of the calculation of the local skin friction as well as the other 
viscous parameters. 

Transition region effect . - The boundary layer method has the capability of 
calculating transition from laminar flow to turbulent flow in two different 
ways. The first approach is to use the transition point as a switching point 
between laminar and turbulent boundary layer calculations. At the transition 
point the turbulent boundary layer calculations are started by activating 
the eddy viscosity coefficient. In general, especially at low Reynolds 
numbers this approach can lead to errors as shown by Cebeci in Reference 28. 
The second approach which is available uses the intermlttancy factor given by 
Chen and Thyson (Reference 29) to modify the eddy viscosity equations to 
account for a region of transition. This modification was developed from the 
point of view of intermittent production of turbulent spots and is a further 
extension of Emmons' spot theory (E.eference 30) . The modification to be used 
is to multiply the inner and outer eddy viscosities equations (16) and (27) 
by the following parameter: 
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( 35 ) 


The effect of this transition region correction can be seen in figure 4 
which compares experimental data to theoretical calculations for local skin 
friction with and without the above correction on a two-dimensional ellipse. 

This transitional effect will be assumed to be the same for axlsymmetrlc bodies. 


Boundary layer transition location . - The location of boundary layer transi- 
tion from laminar to turbulent flow can be either input to the boundary layer 
method or calculated internally within the program. The approach used to 
calculate the transition location is one developed for tiro-dimensional flow 
by Michel (Reference 31) and later verified by Smith (Reference 32). This 
method correlates the local iromentum thickness Reynolds number, Rg and the 
local distance Re 3 molds number, R^, as shown in figure 5 which comes from 
Reference 32. The procedure used is to calculate the values of R^ and Rg 
at each station and to compare them to the curve in figure 5. If the value 
of Rg is less than the value of ^6^ than transition has not been reached 
but if the value of R. is greater than Rfi_,_ then transition has occurred. 

u iR. 


The above method was extended to axisymmetric flow by the use of Mangler's 
transformation. The parameters Rg and R^ are calculated by the axis 3 rmmetrlc 
boundary layer routine and they are then transformed to two-dimensional values 
by the following relationships: 


2-D 




AXISYMMETRIC 


^LOCAL 


X 


2-D 


/ m 


(dX) 


AXISYMMETRIC 


(36a) 


(36b) 
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These values of 6„ and X_ are used to determine values of R. and 
R which can be used in conjunction with figure 5. 

X 

A study of transition location calculation for axis3mtmetrlc bodies was 
recently completed by Kaups (Reference 33) . In this study empirical methods 
due to Granville, Hall and Gibbons, and the method of Michel presented above 
were compared to the stability analysis technique of Smith (Reference 32). 

It was determined that for flows where transition occurred in an adverse 
pressure gradient all of the above techniques predicted transition fairly 
accurately. For flows where transition occurred in favorable pressure 
gradients, only the method of Smith (Reference 32) gave satisfactory results 
as rhown in figure 6 which is taken from Reference 33. The method of Smith, 
however, requires extremely lengthy computer calculation times which makes 
it undesirable for the iterative type of calculation presented in this report. 
Therefore, based on the results of Reference 33, the method of predicting 
transition in the present program should not be used for flows with very large 
Reynolds numbers where the transition location might occur in a favorable 
gradient, but rather the transition point should be input to the program. 

Calculation Procedure 

The viscous flow field about an axlsymmetric body is simulated by calcu- 
lating the Invlscld flow about an equivalent "viscous" body which is formed 
by adding the boundary layer displacement thickness to the original body 
surface. This technique of defining the Invlscid body has been used quite 
successfully for two-dimensional flows as shown in Reference 2 and has 
also been used for axlsymmetric flows as presented in Reference 34. This 
equivalent body is formed by combining the previously discussed geometry 
routine, potential flow method, and boundary layer method under control of 
the axlsymmetric design and analysis method computer program known as ADAM. 

Given the desired axlsymmetric configuration and flow conditions, the 
ADAM program utilizes these sections, as shovm in figure 7, in the following 
iterative manner; 

1. Precise geometry definition for input into the potential flow program. 
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2. 


Calculation of the exact nonlinear potential flow for specified geometry 
and flow conditions. 

3. Calculation of the viscous flow characteristics based on the results 
of the potential flow program. 

4. Addition of boundary-layer displacement thickness to the basic geometry 
for each element. 

5. Recalculation of the pressure distribution utilizing the potential 
flow program, based on the redefined geometry. 

6. Recalculation of viscous flow field based on recalculated pressure 
distribution from redefined geometry, if desired. 

7. Possible iteration of the above scheme; the degree to which this is 
required is presented in the subsequent discussion on correlations 
with experimental data. 


The above technique must be modified when the boundary layer separates 
or when the local body radius approaches zero at the trailing edge of the 
body. When the dimension of the local body radius approaches zero at the trail- 
ing edge, the boundary layer equations become invalid since the 1/r term 
in equation (3) approaches infinity. When this occurs, the boundary layer 
results are ignored from this point downstream to the trailing edge. The 
assumption is then made that the boundary layer displacement area at the 
point where 

6* cos a “ r^ (37) 


is defined by 


DAREA •• TT 



+ 6 * 



(38) 


where p refers to the point where equation (37) is first satisfied. This 
displacement area is then considered to remain constant from the point p 
to the trailing edge. The new "viscous" body coordinates in this region are 
then defined by 


y 

new 



+ DAREA I * 


(39) 


The second problem area occurs when the boundary layer separates from 
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the body creating a separation bubble. This bubble must be accounted for in 

the creation of a "viscous" body if the flow about this configuration is to 

be predicted accurately. The simplest technique of modeling this separation 
bubble ie to assume that the flow leaves the surface parallel to the free- 
stream direction, producing a cylindrical wake shape as shown in figure 8. 

This approach, however , gives decelerations in the flow at the junction of the 
body with the cylinder as shown in figure 9, which do not exist in the real 
flow field. To minimize this problem, a circular arc is used to fair the 
body into the separated cylinder. 

This circular arc is defined by passing a circle through the last three 
"viscous" body coordinates defined prior to the separation point. The radius 
of this circle is then used to create a circular arc which is tangent to the 

"viscous" body at the point of separation. The center of this arc is then 

defined according to whether the surface slope of the body at separation is 
positive or negative. 


If the surface slope is positive then the center is taken as the center 
of the circle passed through the three points as defined above. This center 
is defined by 



dx 


Surface slope at separation 


This arc is then used from the noint of separation to either the end 
of the body or to the maximum point on the arc, where dy/dx = 0, as shown 
in figure 10a. If the maximum point of the arc occurs before the trailing 
edge of the body is reached then a cylinder is defined which extends from the 
maximum point of the arc to the trailing edge. 

If the surface slope is negative then the circular arc is defined such 
that the center is located above the body. The center is then defined by 
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(41a) 


X « X + R Isin tan ^ l|^|] 
c sep L ' 3x'J 

y ^ “ ^sep ^ j^cos tan”^ 


(41b) 


This arc is then used from the point of separation to either the end 
of the body or to the minimum point on the arc, where dy/dx = 0, as shown 
in figure 10b. If the minimum point of the arc occurs before the trailing 
edge of the body then a cylinder is defined which extends from the minimum 
point of the arc to the trailing edge of the body. 

The above separated wake model has been derived from intuitive consider- 
ations rather than from first principals. It does, however, provide reason- 
able results, as will be sho\<m in the subsequent discussion. 

The base drag coefficient for blunt axisymmetric bodies is calculated 
using the method of Hoerner, reference 35. This approach is based on the 
a 88 un 5 >tion that the flow field behind a blunt base is basically a jet pump, 
in that, air flowing aroimd the body leaves the trailing edge forming a 
cylindrical jet which attei^ts to pump away the stagnated air in the base 
region. However , since there is no air to replace this stagnated air, the 
pumping mechanism can only reduce the static pressure acting on the base. 

The effectiveness of this jet punq) mechanism is controlled by the boundary 
layer thickness at the base since this region of lower momentum flow acts 
as a buffer between the stagnated air behind the base and the flow in the 
jet. Since the boundary layer thickness is directly related to the skin 
friction on the body, , Hoerner used Cf to correlate with the base 
drag to develop an empirical approach to determine base drag. Figure 11 
shows the correlation obtained by Hoerner for bodies whose base area is the 
same as the maximum area. This curve is represented by 

*"Dbase (42) 

where the coefficients are based on the base area. Thus, once the skin 
friction on the forebody has been calculated in the boundary layer programs, 
then the base drag can be determined by equation 42. 
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This equation must be modified for boat-tailed bodies, that is, bodies 
whose base area is less than their maximum area. The mechanics of the base 
drag for these configurations do not change, but the calculation must take 
into account the reduced base area. This effect is taken into account by the 
following relationship : 


= C 




and 


so 


BASE 




D 


BASE 


(BOAT TAIL) 


BASE 

BOAT TAIL 




.029 



f 


^base V 

^MAX / 


(42a) 


(43b) 


(43c) 


A comparison of results calculated by the above method in ADAM with experi- 
mental force data from reference 35 is presented In figure 12. One of these 
cases is for a boat-tailed body and the other for a body whose base area is 
also the maximum area. 


The experimental data used for this comparison as well as the configu- 
ration used for the analytical calculations are both subject to some discus- 
sion. The experimental base drag, taken from Figure 4 of Reference 35, 
originally came from an old German report which is not readily available. 

These base drag values were obtained from both force measurements and pressure 
measurements which unfortunately do not agree. Therefore, since it was felt 
that the force measurements were the more accurate, they were used in the 
comparison shown in Figure 12. In addition, no good definition of the con- 
figuration tested was available, therefore, the geometry used in the ADAM 
analysis was taken from the schematics shown in Reference 35. In light of 
these uncertainties the comparison presented in Figure 12 is fairly good in 
that even though the levels are different, the trends are the same. It 
should be noted that this comparison was used only because there is a singular 
lack of experimental data for blunt based axisymmetrlc bodies at low subsonic 
Mach numbers. 
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MPmiMENTAL GOEIELATIONS 


Experimental results from three different configurations were selected 
to establish the extent of validity of the method presented in this report . 

These geometries consisted of a high fineness ratio body of revolution, and 
a sphere in both subcrltlcal and supercritical flow regimes . These corre- 
lations, while limited to some extent by the scope of the present effort, do 
represent a wide range of axlsymmetrlc flow conditions. 

The body of revolution chosen was tested in the low speed wind tunnel at 
the Douglas Aircraft Company, (Reference 36), and is shown in figure 13. This 
model was composed of three sections; an elliptical nose section, a cylindri- 
cal control section, and a parabolic afterbody. The calculation done for this 
configuration used the wind tunnel flow properties, namely, = 71.628 M/Sec 
(235 Ft/Sec), T^ - 288. 3*K (519. 0“R) and Rj^ - 10.05 x 10®. Boundary layer 
transition was fixed on the model and in the calculation at .03048 meters 
(1.2 Inches) from the nose. This model was relatively large for the wind 
tunnel in which it was tested; wall effects, not accounted for in the original 
data reduction, were present. To correct for this, the model was run in the 
potential flow program in the presence of the wind tunnel walls as shown in 
figure 14. The effect of including the walls in the calculation is shown in 
the Inviscld pressure distributions of figure 15. The final results for this 
configuration are shown in figure 16 where the calculated "viscous" results 
are compared to experimental data. The invlscid distribution is also shown 
for reference. In this particular case no separation occurred and so only one 
iteration, that is, two potential flow solutions and two boundary layer solu- 
tions, was necessary. The calculated "viscous" results agree very well with 
the experimental values except in the region of the nose. This discrepancy is 
not due to the calculation method, but rather is due to the model being too 
long for the wind tunnel test section resulting in the nose being in a different 
static pressure field than the rest of the body. The overall effect of viscosity 
on this configuration is seen to be small except in the region of the trailing 
edge. The body is so slender in this region that the boundary layer equations 
are no longer valid so the technique described in the calculation procedure 
was used to modify the viscous body. The results show a pressure oscillation 
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in this wdified region which is due to an «ns»oth curvature distribution. 
However, the level of these pressures agree quite well with the e^erlmental 
values. 

The second case considered was ^ that of a sphere in the supercritical 
flow regime, l.e., ■ 1 x 10®, Since the boimdary layer transition was 

forced to occur at an X/D ■ .65, there were regions of both laminar and 
turbulent flow present. The e]q>erlmental data for this case were taken from 
references 37 and 38. The freestream velocity assumed for this case was 
47.85 M/Sec (157 Ft/Sec). Figure 17 shows the sphere with the "viscous” 
body superln^osed and figure 18 presents a co^arlson between the calculated 
"viscous" solution and eiqierlmental data. Note that while the calculated 
pressure distribution Is in reasonably good agreement with the experimental 
values, the calculated separation point Is ,07 diameters further downstream 
than the e3q>erlmentally measured value. The Invlscld and "viscous" solutions 
for the local skin friction coefficient, G^, are presented in figure 19. 

The "viscious" solution shown Is the fourth iteration, l.e. , the fifth 
potential flow solution, and appears to be the best solution possible for 
this configuration with the technique being used in the present method to 
simulate flow separation. 

The last correlation to be presented is for the flow about a sphere in 
the subcritical regime, l.e. , 8^^ “ 1 x 10®, which is a purely laminar case. 

The experimental data is again taken from Reference 37. The freesteam velo- 
city for this case was assumed to be 4.785 M/Sec (15.7 Ft/Sec). The calcu- 
lated "viscous" body is shown in figure 20 while a comparison of the "viscous" 
pressure distribution to experimental data is shown in figure 21. The calcu- 
lated "viscous" pressures are in close agreement with the experimental values 
with some slight over-prediction in the separated region. The calculated 
separation point is only .03 diameters further downstream than the experi- 
mental value which is excellent considering the large effect that viscosity 
has on this configuration. Figure 22 presents the invlscld and "viscous" 
solutions for the local skin friction coefficient for this case. 
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CONCLUDING lEMAlKS 


A method has been presented for the computation of the viscous flow 
field about axisymmetric bodies at zero angle of attack in incon^ressible 
flow. This con 5 >uting program requires only the specified body geometry and 
desired flow conditions as input . The appropriate theory has been discussed 
and correlations between theoretical and experimental results presented. 

The flow field about axisymmetric bodies at zero angle of attack with 
no flow separation is well defined and can be computed accurately by the 
present method. When flow separation occurs , the flow field Is no longer 
amenable to analytical treatment. Currently , methods do not exist to 
calculate the flow field within a separated region; it is therefore necessary 
to resort to enq>lrical methods to account for flow separation. Since there 
is almost a conq>lete lack of experimental data concerning the behavior of 
separated regions, any en 5 >irlcal methods must necessarily be somewhat crude. 

The most sophisticated model for separation currently available is due to 
Jacob (References 39 and 40) and is strictly for two-dimensional airfoils. 

An unsuccessful attempt was made In Reference 41 to adapt Jacob's approach to 
axisymmetric configurations. The conclusions of Reference 41 indicated that 
the assumed boundary conditions needed to be modified if this approach was to 
be used for axisymmetric flow. It is proposed that the Douglas-Neumann 
program be used to pursue this approach at modeling separation. This 
potential flow program is ideal for attempting to use Jacob's technique since 
it already has the ability to specify a non-uniform flow distribution over 
all or part of a configuration; therefore, only suitable boundary conditions 
would have to be added to the program. It is felt that this approach can be 
successful in modeling separation if care is taken in developing the distri- 
bution of non-uniform velocity as well as specifying the proper boundary 
considerations . 

The further extension of this model to the calculation of flow about 
axisymmetric bodies at angle of attack is also possible. The potential flow 
routine contained in the present method has the capability of predicting the 
flow field about non- Ilf ting bodies at angle of attack by combining the stream- 
flow and the crossflow solutions. The boundary layer analysis would require 
the replacement of the routine in the present method by a three-dimensional 
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technique, which is currently not available . However , it is felt that a 
good approximation to the boundary layer calculations can be made by the 
small crossflow program of Reference 42 , 

One area of primary concern in extending the method to include an angle 
of attack capability is the determination of the separation line about the 
body. The present method of predicting separation for two-dimensional bodies 
and for axisymmetrlc bodies at zero angle of attack is to find the location 
where the skin friction goes to zero. It has been shown in several studies, 
Including those reported In References 43 and 44, that this condition does not 
apply in three-dimensional flows because the skin friction along a separation 
line is not necessarily zero. Therefore, some method of determining the 
separation line for axisymmetrlc bodies at angle of attack must be developed. 
It is proposed that the present method could be extended to calculate the 
"viscous" flow about axisymmetrlc bodies at angle of attack when no flow 
separation is present. This method could then be used to assist in the 
development of a procedure for determining the separation line location. Once 
the location of the separation line is known then a model could be developed 
for analyzing the viscous flow about the separated body. The development of 
such procedures is not a simple task and considerable effort would have to 
be expended; but the reward for accon?>lishing this task is an advance in the 
ability to calculate the real flow about arbitrary three-dimensional bodies 
which is our ultimate goal. 
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APPENDIX A 


INPUT INFORMATION FOR ADAM COMPUTER PROGRAM 

This part of the report contains the necessary information to input 
data to the ADAM computer program. The input data is broken into three 
sections: smoothing, potential flow, and viscous flow. These sections can 

be used together in the iterative fashion described in the main text, or the 
potential flow and viscous Flow sections may Ije used independently. A 
detailed card-by-card description of all input quantities is given follov;ed 
by a set of input forms which can be. used to facilitate the loading of tlie 
input data into the program. 


26 



Input Instructions 


The Adam program requires one system control card followed by the 
required sets of data cards for each program option to be executed. The 
sets of data furnished must be In the same order as the options are specified 
on the system control card. If an Iteration Is desired the system control 
card Is repeated along with the necessary other data cards. 

The general scheme used In describing the Input data Is shown below: 

Routine 

Column Code Format Explanation 

Column - Column Indicates the starting position on the card for each 
data field. 

Code - The "code" gives the FORTRAN name used In the read statement by 
the program. 

Routine - "Routine" Indicates the subroutine where the data Is read. 

Format - The parameter "FORMAT" which is given right under the routine 
name, indicates the FORTRAN format of the data read statement 
field. The parameter 15 would indicate that the parameter is 
an Integer In a field that is 5 columns wide. Integers should 
be punched on the right side of the field (right justified). 

The parameter FIO.O would indicate a fixed point number punched 
with a decimal point (i.e., -12.354). The number may be punched 
anjn/here in the field Indicated irrespective of the decimal point 
location Indlcatied by the format. The parameter E12.6 would 
indicate a floating point number punched with a decimal point 
(i.e,, 5.0 X 10^). The number must be punched to the right 
of the field in the manner 5.0E+06. 

Explanation - The description of the input data is given under "explanation". 
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SUSTEH CONTROL DATA CARD (This card must be the first card in the data deck) 


Column 


4 


8 


12 


16 


17-20 


Code 

IGEOM 

INEUM 


IBOUND 


ITER 

IFINSH 


Routine 

Format 

MAIN 

II 


MAIN 

II 


MAUI 

II 


MAIN 

II 

MAIN 

14 


Explanation 
Smoothing option flag 

■0 no smoothing is desired 
“I smoothing is desired 
Potential flow option flag 

■0 No potential flow solution is 
desired 

■1 Potential flow solution is 
desired 

Boundary layer option flag 


■O No boundary layer solution is 
desired 

■1 Boundary layer solution is 
desired 

"Viscous" body formation flag 

■•0 No "viscous" body is formed 
-1 "Viscous" body is formed 
Termination Flag 

■0 Another case expected 

■9999 Program will stop after exer- 
cising all options specified 
above 
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SMOOTHING SECTION 


These cards required 

if IGEOM 

» 1 on system control card. This section is 

used to smooth body geometry data before it is input to the potential flow 

program. 




Smoothing Control Card 



Column 

Code 

Routine 

Format 

SMOOTH 

Explanation 

2 

NPTS 

Number of input data points for this 



13 

configuration. 

NPTS must be < 100 

8 

ITAPE 

SMOOTH 

Data source flag 



11 

=0 Data input on unit 5 (card input) 
1*0 Data input on unit 1. This is 




used for a case where a "viscous" 
body generated by the iteration 
procedure is being read. 

Geometry Datd Input Cards 



These i 

cards are input only if 

ITAPE “ 0. 

Column 

Code 

Routine 

Format 

Explanation 

x-coordinate 

cards 



1-10 

x(l) 

SMOOTH 

x-coordlnates starting at the leading 

11-20 

x(2) 

6F10.0 

edge and proceeding along the upper 

21-30 

etc. 

x(3) 


surface to the trailing edge. Input 
6 x-values on each card. The numbers 
of x-values must be equal to NPTS. 

y-coordinate 

cards 



1-10 

y(i) 

SMOOTH 

y-coordlnates to correspond to the 

11-20 

y(2) 

6F10.0 

above x- locations, y-values must 

21-30 

y(3) 


be positive. Input 6 values per 

etc. 



card. 
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POTENTIAL FLOW SECTION 


These cards required if INEUM ■= 1 on system control card. The Input 
geometry for this program may be obtained from the geometry storage unit (10) 
as generated by the smoothing section, or It may be Input directly on unit 5. 
Thus, this program may be operated as a separate entry If so desired. The 
program saves the geometry data element midpoints with the corresponding 
pressure coefficients on unit 3 for input to the boundary layer routine and 
it saves the basic non-dimentional input Neumann coordinates on unit 1 for 
use if a "visous ' body is desired. 

Title Card 


Column 

Code 

Routine 

Format 

Explanation 

1 

HEDR 

PARTI 

10A6 

Title of case. May be any characters 
input in the first 60 columns of card 

63 

CASE 

PARTI 

16 

Case number 

77 

PSF 

PARTI 

16 

Additional identifier for this case. 

Flag Card 




Card 

columns 

1-30 when punched with any non-zero integer, activate 

flags that 

indicate the following: 


Column 

Code 

Routine 

Format 

PARTI 

11 

Explanation 

1 

NB 

The number of bodies input. Normally 
set equal to 1. 1 < NB £ 5 

2 

NNU 

PARTI 

11 

The number of non-uniform onset 
flows. Normally set equal to 0. 

3 

FLG03 

PARTI 

Axisymmetric flow flag. 



11 

■0 No axisymmetric stream-flow 
solution calculated. 




■1 Axisymmetric streamflow solution 
is calculated 




Normally set equal to 1 
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Flag Card 
Colunm 
4 


5 


6 


7 

8 


11 

12 


(Continued) 

Code 


Routine 

Format 


FLG04 PARTI 

II 


FLG05 PARTI 

II 


FLG06 PARTI 

II 


Explanation 

Cross flow flag. 

=0 No cross flow solution is 
calculated 

-1 Cross flow solution is calculated 

Normally set equal to 0 

Off-body point flag 

=0 No off body points input 

"1 Off body points are input 

This flag allows the velocity at 
points off the body surface to be 
determined. 

Basic data formation flag 
=0 A full case will be done 

=1 The basic data, i.e., midpoints, 
normals, etc. will be formed and 
printed. No velocities will be 
calculated. 


FLG07 PARTI 

II 


FLG08 PARTI 

II 


FLGll PARTI 

II 


FLG12 PARTI 

II 


Ellipse generator option 

=0 Body coordinates will be input 

=1 An ellipse is generated using data 
input later. No body coordinates 
are input 

Matrix print flag 

=0 Coefficient matrices are not 
printed. 

"«1 Coefficient matrices will be 
printed. 

Normally set equal to 0 

Perturbation velocity flag 

=0 Normal case 

=1 No onset flow used. Only per- 
turbation velocities are calculated. 

Potential matrix solution 

=0 Normal case 

=1 A potential matrix is solved 
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Flag Card (Continued) 


Column 

Code 

Routine 

Format 

Explanation 

13 

FLG13 

PARTI 

11 

Matrix solution flag 
“0 No matrix solution done 
“1 Matrix solution performed 
Normally set equal to 1. 

14 

FLG14 

PARTI 

11 

Prescribed tangential velocity flag * 
■0 Normal case 

=1 Tangential velocities are specified 

15 

FLG15 

PARTI 

11 

Strip ring vortlcity flag * 

■0 Itormal case 

■1 A vortlcity distribution is 
formulated. 

16 

FLG16 

PARTI 

11 

Axisymmetric uniform flow flag 
«0 Normal case 

“1 Axisymmetric uniform flow solution 
is omitted 

Normally set equal to 0. 

17 

FLG17 

PARTI 

11 

Crossflow uniform flow flag 
“0 Normal case 

■1 Crossflow uniform flow solution 
is omitted. 

Since FLG04 is normally 0 then 
so is FLG17 normally set equal to 0. 

18 

FLG18 

PARTI 

11 

Surface vortlcity flag * 

■bO Normal case 

“1 Surface vortlcity is generated. 

19 

FLG19 

PARTI 

11 

Prescribed vortlcity Flag 
B>0 Normal case 

■1 A prescribed vortlcity is input 

20 

FLG20 

PARTI 

11 

Total vortlcity flag 
■0 Normal case 

■1 Total vortlcity calculated 



Flag Card (Continued) 


Column Code 


Routine 

Format 


Explanation 


21 


22 


23 


28-29 


FLG21 


FLG22 


FLG23 


NIN 


PARTI 

11 


PARTI 

II 


PARTI 

II 


PARTI 

12 


Extra crossflow flag * 

=0 Normal case 

=1 Extra crossflow option used 
Generated boundary condition flag * 

=0 Normal case 

“I Boundary conditions generated 
Ring wing option flag * 

=0 Normal case 

=1 Ring wing option used 

Tape input flag 

=0, 10 Data input on unit 10 from 
smoothing program 

=5 Data input from unit 5 (card input) 


30 ITER PARTI Iteration tape flag 

II 

=0 x/c, y/c transformed data saved 
on unit 15 

=1 x/c, y/c transformed data not 
saved . 

This flag is necessary because for 
a "viscous" body to be formed, the 
coordinates of the original unmodified 
body must be saved. Therefore, for 
the first case set ITER = 0. For 
subsequent Iterations we do not v/ant 
to use the modified bodies to form 
new bodies so set ITER = 1. 


* These flags are for special options x^hlch are discussed in the main 
text. They are never used for a normal axisymmetric calculation. 
Therefore, set them equal to zero. 


33 



Chord Card 


Column 

Code 

Routine 

Format 

PARTI 

FIO.O 

Explanation 

1 

CHORD 

Reference chord length used to 
non-d imens lonal 1 z e x and y 
coordinates 

11 

MN 

PARTI 

FIO.O 

Mach number (MN < 1.0) use to 
approximate effect of compressi- 
bility (Gothert's rule) 

21 

TCNST 

PARTI 

FIO.O 

This Is a cons tant which is used 
for the value of the tangential 
velocity If this option la desired. 

Body Transformation Card 



Column 

Code 

Routine 

Format 

BASICl 

13 

Explanation 

8 

NN 

The number of Input points on this 
body. NN S 100 

11 

MX 

BASICl 

FIO.O 

A factor used to multiply all 
x-coordlnates. MX Is assumed equal 
to 1 If no value Is Input. 

21 

MY 

BASICl 

FIO.O 

A factor used to multiply all y- 
coordlnates. MY Is assumed equal 
to 1 If no value la Input. 

31 

THETA 

BASICl 

FIO.O 

An angle (in degrees) through which 
all points of a body are to be 
rotated about the origin in the 
clockwise direction. 

41 

ADDX 

BASICl 

FIO.O 

A constant to be added to all 
x-coordlnates 

51 

ADDY 

BASICl 

FIO.O 

A constant to be added to all 
y-coordinates 
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Body Control Card 


Colunm 

Code 

Routine 

Expl 

anation 







10 

BDN 

BASICl 

11 

Body sequence number. Tills program 
will handle up to 5 bodies. 

20 

SUBKS 

BASICl 

11 

Subcase Flag. 

=0 Normal case 

“1 Use unmodified coordinates of 
the previous case. 

30 

NLF 

BASICl 

11 

Non-lifting flag 

=0 Body is non-lifting 
(normal case) 

•=1 Body is lifting (this is used 
in special’ option) 

31 

XE 

BASICl 

FIO.O 

Value of major semi-axis for use 
by ellipse generation option. 

41 

YE 

BASICl 

FIO.O 

Value of minor semi-axis for use by 
ellipse generation option 

Note: if XE = YE a sphere will 

be formed. 


Geometry Data Cards 

The body geometry data cards are included only if the input parameters 
NIN “ 5 and FLG07 “ 0 on the flag card. If NIN =0 or 10 then the 
data is read from unit 10. If NIN = 5 and BDN = 0, then the following 
cards contain the x-y coordinates of off-body points instead of x-y 
geometry data. The number of either geometry data point or off-body points 
must be equal to NN. 

x-Coordinate cards (six values per card) 


Column 

Code 

Routine 

Format 

Explanation 

1 

TXl(l) 

BASICl 

x-coordlnates of body input from 



6F10.0 

leading to trailing edge. 

11 

TX1(2) 



21 

TX1(3) 
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y-Coordinate 

cards 

(six values per card) 


Column 

Code 

Routine 

Format 

Explanation 

1 

TYl(l) 

MSICl 

6F10.0 

y-Coordinates of body which 
correspond to the x-values ; 
y values must be positive. 

11 

TY1(2) 



21 

TY1(3) 



etc. 





NOTE: Each body input, including the off body points, requires the body 

transformation card, the body control card, and may also require the 
geometry data cards depending on the input flags. This is the 
stopping place for a normal axisjrmraetric case. The following cards 
are input only if one of the special options is required. 


Tangential Velocity Data (six values per card) 

These cards are input only if FLG14 ^ 0 and TCNST =0.0 
Column Code 

1 TG(1) 

11 TG(2) 

21 TG(3) 

etc. 

Non-uniform Flow Cards (six values per card) 

These cards are input only if NNU ^ 0. 


Column 

Code 

Routine 

Format 

Explanation 

6 

NUM 

BASIC2 

15 

Non-uniform flow identification 
number. 


Routine 

Format 


Explanation 


BASICl Specified tangential velocities at 

6F10.0 element midpoints. 
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Non-uniform Flow Cards (Continued) 


Column Code 


Routine 

Format 


19 MSF BASIC2 

12 


21 TYPE BASIC2 

FIO.O 


31 FG BASIC2 

FIO.O 


Explanation 


If MSF » 0 the flow velocities 

Nq,Tq will be used for the axisym- 
metrlc case only. 

If MSF “ 1 the flow velocities 
Nq.Tq will be used for the cross 
flow case only. 

If MSF > 1 the flow velocities 
will be used for both axisymmetric 
and cross flow cases. 

Flag which specifies the type of input 
flow velocities at each mid-point. 

If TYPE > 0.0, the velocities are 
input as x & y components. 

If TYPE = 0.0 the velocities are 
input as normal & tangential 
components. 

If TYPE < 0.0 the automatic generation 

of the flow due to a rotating body 
is used. 

Constant used by the flow generator. 
Type must be less than 0.0. 


The following cards are input only if NNU ^ 0 and TYPE 1.0. 
Normal velocity cards (six values per card) 


Column 


Code 


N0(1) 


Routine 

Format 

BASIC2 


Explanation 

This is either the x or normal 
velocity component depending on the 
value of type above. These values 
must be in sequence with the coordi- 
nate data. If the x component is 
input it is defined as positive to 
the right. If the normal velocity 
is input it is positive if it is 
to the interior of the body. 

NN-1 values are input. 
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(six values per card) 


Column 


Code 


Routine 

Format 



1 T0(1) BASIC2 This Is either the y or tangential 

6F10.0 velocity component depending on the 

value of type above. These values 
must correspon to the NO values above. 
If the y component is input it is 
defined as positive if it is orientated 
upwards . If the tangential velocity 
is input it is positive if the flow 
field is to the left of the vector 
representing the tangential velocity. 
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VISCOUS FLOW SECTION 


These cards required if I BOUND = 1. 


BOUNDARY LAYER- PROGRAM 

The geometry and pressure distribution data required by this program 
may be input directly on cards (Unit 5), or read from the data save unit 
(Unit 3) as generated by the Neumann program. 


HEADER CARD 


This card is supplied purely for description purposes. 
Column Code 


1-60 TITLE 

61 CASE 


Flag Control Card 

This card contains flags which control the type of flow to be considered 
and the form of the input. 

Explanation 


The number of the x-statlon where the 
flow goes turbulent measured from the 
stagnation point (l.e., the leading 
edge for axisymmetric bodies at zero 
angle of attack) if transition is to 
be calculated by the program set NXT 
to be one greater than the number of 
points input. 

Transition flag 

=0 Boundary layer transition point 
is input 

=1 Boundary layer transition point 
is computed. Set NXT to be 
greater than number of points 
input . 


Column 


Code 


Routine 

Format 


NXT 


INPT 

14 


LG16 


INPT 

II 


Routine 

FoTnaat 

Explanation 

INPT 

15A4 

Description of input 

INPT 

Case number 

A4 
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Control Card (Continued) 


Colunm Code 


Routine 

Format 


Explanation 


6 LG17 


7 LG18 


8 LG32 


9 LG26 


10 LG40 


11 LG41 


INPT 

II 


INPT 

II 


INPT 

II 


lOTT 

II 


INPT 

II 


INPT 

II 


Transition control flag 

“0 Transition is instantaneous 

=1 Transition is gradual (transitional 
region used) 

Transverse curvature flag 

=0 No transverse curvature correction 
used. 

=1 Transverse curvature corrections 
applied . 

Print control flag 

=0 Print using long format (with 
velocity profiles) 

=1 Use short printout (no velocity 
profiles) 

Velocity input control flag 

=2 Velocity ratio (Ug/U^) is input 

=3 Pressure coefficient (c ) is 
input . ^ 

Unit input flag for geometry and 

velocity data. 

-0 Data read from unit 3 as 
generated by the potential 
flow program. 

/O Data read from cards (unit 5) 

System of units FLAG 

=0 English system of units 

=1 Internation system of units 


Flow Condition Card 


Column Code 


Routine 

Format 


Explanation 


1 TI 


INPT Reference static temperature used to 

FIO.O compute the reference fluid proper- 

ties. If TI is input as zero then 
TI is set equal to either 288.33°K 
or 519° R depending on FLAG LG41 
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Flow Condition Card 

(Continued) 


Column 

Code 

Routine 

Format 

INPT 

FIO.O 

Explanation 

11 

RMI 

Reference or free-stream Mach 
number. 




-0.0 UI is input next 

/0. 0 UI is computed from PvMI. 

21 

UI 

INPT 

FIO.O 

Reference or free-stream velocity 

=0.0 is input above 

7 * 0.0 input as zero above. 

31 

FK 

INPT 

FIO.O 

Flow index 

=0.0 2-D flow assumed 

-1.0 Axlsyinmetrlc flow assumed 

41 

RL 

INPT 

FIO.O 

Chord or reference length 

51 

RI 

INPT 

E12.0 

Reynolds number/foot 

, Uoo 

R /«, = — 

C V 




If CHORD - 1.0 then RI must be 
Reynolds number based on CHORD. 




NOTE! The input of either Mach mmber 
or freestream velocity is for 
convenience only. This program 
is entirely incompressible. 

Radius card 




Column 

Code 

Routine 

Format 

Explanation 

1 

ROMAX 

INPT 

FIO.O 

Maximum radius of body. This is used 
to obtain frontal area for skin 
friction calculation. 

11 

DETAl 

INPT 

FIO.O 

Initial step size of boundary layer 
velocity profile grid. For a case 
which contains turbulent flow set 
DETAl •= .005. 
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Radius card (continued) 


Column Code 


Routine 

Format 


Explanation 


21 VGP 


INPUT 

FIO.O 


VGP is the growth factor for the 
boundary layer velocity profile grid; 
for cases with, turbulent flow set 
equal to 1.14. 

NOTE: For laminar cases the boundary 

layer velocity profile grid 
may be made constant if VGP = 1.0 
is input. However, if this is 
done the minimum value of DETAl 
that can be input is approxi- 
mately .10. This can be calcu- 
lated if the value of the 
transformed boundary layer 
thickness, ETAINF, in known. 

Then DETAl becomes 


DETAl 


ETAINF 

100 


Geometry-Pressure Distribution Cards 

These cards input only if LG40 0. 


Point Number 

Card 



Column 

Code 

Rout ine 
Format 

Explanation 

1 

NXM 

INPT 

14 

Number of data points to be input. 
Maximum of 100 points allowed. 

x-Coordinate 

Data Cards 



Column 

Code 

Routine 

Format 

Explanation 

1 

11 

21 

etc. 

y-Coordinate 

XS(1) 

XS(2) 

XS(3) 

Data Cards 

INPT 

6F10.0 

x-coordinate points input from 
leading to trailing edge input 6 
points per card. Number of points 
=■ NXM 

Column 

Code 

Routine 

Format 

Explanation 

1 

YS(1) 

INPT 

6F10.0 

y-coordinate points corresponding 
x-coordinates above input 6 points 
per card. 
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y-Coordinate Data Cards (continued) 


Column Code 


Routine 

Format 


11 YS(2) 

21 YS(3) 

etc. 


Pressure Distribution Cards 


Column Code 


Routine 

Format 


I JKS.o 

II UE(2) 

21 UE(3) 

etc. 


Explanation 


Explanation 


Velocity-pressure-distribution points 
corresponding to x-polnts Input above 
input 6 points per card. 

If LG26 - 2 Ue/U^ input 
LG26 - 3 Cp input 
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AXISYM^^ETRIC DESIGN AND ANALYSIS METHOD 
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Instructions to Keypunch: ENGINEER PHONE 

Do not punch blank columns 

DATE PAGE 



SMOOTHING PROGRAM 
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HEADER CARD NEUMANN POTENTIAL FLOW PROGRAM 
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ENGINEER PHONE 

Instructions to Keypunch 

Do not punch blank columns DATE PAGE 
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THIS INPUT FORM IS THE SAME FOR GEOMETRY DATA INPUT TO THE SMOOTHING ROUTINE, 
THE POTENTIAL FLOW PROGRAM, AND THE BOUNDARY LAYER ROUTINE. THIS FORM IS ALSO 
USED FOR VELOCITY DATA INPUT TO THE BOUNDARY LAYER PROGRAM AND FOR THE NON- 
UNIFORM VELOCITY DATA WHICH CAN BE INPUT TO THE POTENTIAL FLOW ROUTINE. 




APPENDIX B 


CONTROL INFORMATION FOR ADAM COMPUTER PROGRAM 

This part of the report contains the necessary control information to 
operate this program on a computer system. This section contains the overlay 
structure as well as flow charts of the main subroutines including input flow 
information. Also, the various data sets used between main programs are 
described. 
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OVERLAY STRUCTURE OF ADAM 


MAIN 

NEUMAN 

INSl 

ARSIN 


PARTI 

BASICl 

BASIC2 

MATRIX 

NOTS 

XYZ 

QC 

HLAMB 

XYZl 

XYZ2 

ELLC 

INEL 

EL IP 

ELINT3 


PREP 

SOLVIT 



OVERLAY C0,Q) 

OVERLAY 

Cl,Q) 

OVERLAY 

(3,0) 

OVERLAY 

,C4,0) 

OVERLAY 

(5,0) 

OVERLAY 

(6,0) 


PART4 


OVERLAY 

C4,l) 


BOUNDL 

EINF 

HEAD 

IVPF 

FLPR 

EDVS 

SHFT 

MOMC 

OTPT 

TENS 

INPT 

SLOPE 


SMOOTH 

SM5PT 


OVERLAY 

(7,0) 


ITERAT 

TABLE! 

CIRCLE 


OVERLAY 

(4,2) 


OVERLAY 

(4,3) 


AXIS 

COMBO 

SOLCOM 


CROSS 


EXCROS 
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BAglG FLOW CHART FOR SUBROUTINE PARTI 



5 ' 










BASIC FLOW CHART FOR SUBROUTINE BASICl 


START 


READ BODY TRANSFORMATION CARD FROM UNIT 5 


© 


READ BODY CONTROL CARD FROM UNIT 5 



/'1F\ 

PLG07 


READ BODY COORDINATES 
FROM UNIT 13 

(SAVED FROM PREVIOUS CASE) 
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BASIC FLOW CHART FOR SUBIOOTIIE BASICl (CONTINUED) 










































BASIC FLOW CHART FOR SUBROUTINE BOUNDL (CONTINUED) 
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BASIC FLOW CHART FOR SUBIOHTIIE BOTODL (CONTINUED) 



yes 


CALL SHFT 
CALL FLPR 
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The following is a description of the output symbols from the various 
sections of the ADAM computer program; 


SYMBOL 

ADDED MASS 
AJK 


ADDX 

ADDY 

BJK 

BODIES 
BODY NO. 

CHORD 

COSA 

CP 

DALPHA 

DELTAS 
MACH NO. 
MX 

MY 


NEUMANN POTENTIAL FLOW SUBPROGRAM 

DESCRIPTION 

E2tt . Phi • Vj^.ds 

Influence coefficients resolved parallel to the outward 

normal of the element . Output only if FLG08 «• 1. 

Constant which is added to all X-coordinates of a particular 
body. Value printed out for each body. 

Constant which is added to all Y-coordinates of a particular 
body. Value printed out for each body. 

Influence coefficient resolved parallel to the tangent 

direction of the element. 

Number of bodies in system, same as NB input on flag card. 

Number of this particular body. This parameter input on body 
control card. 

The reference chord for the system. 

The cosine of DALPHA 

The pressure coefficient on a body element. 

The change in angle between consecutive elements of a body, 
(degrees) 

The length of a body element. 

Mach number used in Gothert’s transformation. 

The factor by which all X-coordinates are multiplied for one 
body. Input on body transformation card. 

The factor by which all Y-coordinates are multiplied for one 
body. Input on body transformation card. 
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SYMBOL 

N 

NN 

NNU 

PHI 

PSF NO- 
SIGMA 
SINA 
SUM(T) 

SUMDS 

TCNST 

THETA 

TI 

VOLUME 

X 

XE 

Y 

YE 


DESCRIPTION 

Velocity normal to a surface element, this is a measure of 
how well the boundary condition of zero normal velocity is 
satisfied. 

The number of geometry data points for a given body. This is 
input on the body transformation card. 

The number of non-uniform onset flows to be considered. 

Value of potential on each surface element. 

Identification for this case. . 

Source density on each surface element. 

Sin of DALPHA 

DELTA(S) This is the sommation of T multiplied by Deltas up to each 
element midpoint. 

Summation of Deltas, surface distance around the body. 

Constant value of tangential velocity used in special option. 

The angle through which a body is to be rotated about the 
origin in a clockwise direction. 

The velocity at each midpoint. 

The volume of the body being analyzed, (calculated by Neumann) 

The input X-coordlnate defining the body surface, or off-body 
X-coordinates. 

The value of semi-major axis used in ellipse generation option. 

The input Y-coordinates defining the body surface, or off-body 
Y-coordinates. 

The value of semi-minor axis used in ellipse generation option. 
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OUTPUT DATA SYMBOLS 


Finite Difference Boundary Layer Subprogram 

Output of this routine consists of CASE DATA and STATION DATA Inputs as 
well as the computed STATION DATA. Body geometry data, flags and counters , 
and reference quantities are printed out under the heading of CASE DATA. 
Values of parameters at the outer edge of the boundary layer as well as the 
boundary condition inputs are printed out under the heading of STATION DATA. 
These are followed by iteration results, velocity profiles for each x-station 
(if FLG32 = 0), and a summary of the computed boundary-layer parameters as 
functions of streamwlse or x-dlstance. 

Error messages generated by the program are printed out at the end of 
the STATION DATA printout if they are generated by input errors. Other error 
messages are Issued at different locations in the profile printout if errors 
are detected during the computations. 


SYMBOL 

ALPHl 

ALPH2 

BETA 

C 

CDBASE 

CF 

CFA 

S 

DELS 

DELVW 

EPS 

EPSl 

ETA 


DESCRIPTION 

Local body slope dy/dx. 

Not Used in this program. 

B = (25/u^)(du^/dO 
CHORD 

Base drag coefficient. 

= t^/(1/2 value of local skin friction coefficient. 

Total integrated skin friction to each point. 

Pressure coefficient 

Boundary layer displacement thickness. 

Delta V(l,2) used in iteration for v(l»2) . 

E , eddy viscosity parameter for outer region. 

Ej^, eddy viscosity parameter for inner region. 

q, non-dimensionalized boundary layer thickness to each 
point in the boundary layer. 
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SYMBOL 


DESCRIPTION 


CTAE 

ETMKF 

F,FPjFPP 

FPPW 

FPW 

FW 


GW 

GPW 

H 

HE 

HI 

IMAX 

K 

KK 

ME 

MUE 

MREF 

MUREF 

PE 

PRO 

QW 

REY 


value of n which corresponds to 6. 

Non-dimensional boundary layer thickness used as maximum, 
value in forming numerical solution grid. 

f ,f ’,f", respectively, 
f" at the wall, 
f * ■ U/U at the wall. 

f , this is the transformed stream function at the wall, 
w 



Not used in this method. 

Not used in this method. 

Boundary layer form factor, H » 6*/6 
Enthalpy 

Initial step size, same as DETAl in input. 

Number of points taken through the boundary layer. 

Initial step size of the variable grid system. 

Variable grid parameter chosen internally. 

Local Mach number. 

Local dynamic viscosity, at edge of boundary layer. 

Free stream Mach number. 

Free stream dynamic viscosity, 

Pressure at edge of boundary layer, P^. 

Laminar Prandtl number. 

Not used in this program. 

Reynolds number based on reference conditions (see input) 
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SYMBOL 


DESCRIPTION 


RHOREF 

R„/C 

RR 

RTHETA 

S 

S/C 

SHORT? 

SQUIG 


ST 

SWEEP 

TE 

THETA 

TREF 

TRFLAG 

TRINT 

TW 

TVC 

UE 

UPLUS 


Freestream reference density. 

Rq/L axlsynmetric radius. 

Not used by this program. 

Momentum thickness Reynolds number, Rg. 
Reynolds number based on local conditions. 



Surface distance. 

Nondlmenslonallzed surface distance. 

Flag which tells program to print velocity profiles. 
Same as FLG32 in input. 

Transfoirmed x-coordinate, 5 



Not used in this program. 

Not used in this program. 

Temperature through boundary layer. Not needed for this 
program. 

Momentum thickness, 6. 

Reference temperature, T^. 

Flag which determines transition (input). 

Flag which determines instantaneous transition or use of 
transitional region option (input). 

Temperature at the wall. Not used in this program. 

Transverse curvature flag (input). 

Velocity at edge of boundary layer. 

Non-dlmenslonlized velocity in the boundary layer. 
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SYMBOL 


DESCRIPTION 


X 

X/C 

XI 
Y 

Y/C 

YPLUS 

VREF 


X-coordinate 

Non-dimens ionalized x-coordinate 

Transformed x-coordinate - same as SQUIG 
Y-coordinate 

Non-dimenslonalized y-coordinate 

Non-dimenslonalized y-coordinate in boundary layer. 
Reference velocity (input) 


Iteration Subprogram 

XNEW and YNEW These are the coordinates of the equivalent viscous 
body. The original coordinates modified by the addition of the boundary 
layer displacement thickness 6 *. 
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DESCRIPTION OF STORAGE UNITS 


Hie 

TAPE 

TAPEl 

TAPE2 

TAPES 

TAPE4 

TAPES 

TAPE6 

TAPES 

TAPE9 

TAPEIO 


following Is a description of all disk storage units used in ADAM; 

DESCRIPTION 

This unit used in subroutine SOLVIT as a scratch unit 
and also used to transfer the "viscous" coordinates 
from subroutine iterat to subroutine smooth or sub- 
routine BASICl, 

This unit used in subroutine SOLVIT as a scratch unit 
and also used to transfer the boundary layer displace- 
ment thickness's from subroutine OTPT to subroutine 
ITERAT. 

This unit used to store source densities in subroutine 
SOLVIT and used to transfer body geometry and pressures 
from subroutine AXIS to subroutine INPT. 

sin 01 , cos a, TCNST, TG<I), cos R^, 2| (sin a) (cos a) | , 
“o* "^o' ^N» V ®tc. 

Used exclusively in Neumann subprogram to store and 
transfer data. 

This tape used for card input. 

This tape used for printed output. 

Tills tape used to store extra cross flow matrices, 

EC , ECY, ECZ, in subroutine MATRIX. 

This tape used to store axisymmetrlc flow matrices 
AS, AY, AZ, in subroutine MATRIX. 

This tape used to store cross flow matrices CX, CY, CZ 
in subroutine matrix and also used to transfer smoothed 
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TAPE 


DESCRIPTION 


TAPEIO 

TAPEll 

TAPE12 

TAPE13 

TAPE15 


(Continued) coordinate data from subroutine smooth to subroutine 
BASICl. 

This tape used as a scratch unit in subroutine SOLVIT. 

This tape used to store transformed parameters XI, Yl, 
X2, Y2, and AS^. 

This tape used to store untransformed coordinates 
(TXl, TYl) for use in SUBCASE option. 

This tape used to store transformed coordinates, 

(XI, Yl) for use in subroutine HERAT. 
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APPENDIX C 


PROGRAM LISTINGS 


This part of the report contains the source card listings for the 
axisymmetrlc design and analysis method (ADAM) computer program. This program 
may be run either on a CDC or an IBM computer. The listing as presented here 
is for the CDC version of the program. This program has been run on the CDC 
6600 computer. The program is written in FORTRAN for the CDC run compiler 
and has been run under the scope 3.1 and 3. A operating systems. In this 
listing all cards that are peculiar to the CDC version of FORTRAN are identi- 
fied by a C in card column 80. All cards that are peculiar to the IBM FOR- 
TRAN IV compiler are identified by an I in card column 80 and a C in card 
column 1. In other words, the code for both CDC and IBM machines is in the 
deck but the IBM cards are made inactive by converting them to comment state- 
ments (C in card column 1). Since all of the machine dependent cards are 
Identified by an I or C in card column 80 it is a simple matter to convert 
the deck from one version to the other with a small conversion program. When 
converting from CDC to IBM code this conversion program reads and copies each 
card to a storage unit. If a card has a C in card column 80, then a C 
is written into card column 1 to make the CDC peculiar card inactive. If a 
card has an 1 in card bolumn 80, then the C is removed from card column 1 
and the card image written to the storage unit as an active FORTRAN statement. 
The conversion from IBM back to CDC is made in a similar manner. The con- 
version program to convert the deck from CDC to IBM FORTRAN is listed below 
(for use on an IBM machine) : 

DIMENSION DATA(22) 

DATA CB,CC,CI/1H, 1HC,1HI/ 

REWIND 19 
DO 100 1=1,20000 

READ (5,20,END=300) DATA 
20 FORMAT (1A1,19AA,1A2,1A1) 

IF (DATA(22) .EQ. C7) DATA(l) = CB 
IF (DATA(22) .EQ. CC) DATA(l) = CC 
WRITE (19,20) DATA 
100 CONTINUE 
300 STOP 
END 
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This program places the new deck with IKI cards made active, and CDC cards 
inactive, on to unit 19. The references to unit 19 above can be changed 
to unit 7 to punch the deck out . 
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17 ir? .NOT, RSMALL) GO TO 19 XYZ2 075 
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NSW ■ I XYZ2 101 

C*** ***CALCULATE quantities which have not yet been calculated for Isfl XYZ2 102 
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Ro rn isao comb 1 76 

1750 CONTINUE: COMB 177 

return comb 17 « 

END COMB 179 
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SUBROUTINE SOLCOH( OV , A . ICNj) SOLC 001 

SOLC 002 

note that this subroutine solves UnLy a 1X1 or 2X? matrix soli 005 

A* it is separated so that IE program Is ever INLARGEP, this sotc ooa 

** TS WHERE THE MATRIX SOLUTION EOR ThE COMBINATION PART OF SOLC OOS 
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CROS 032 

dimension VX(lOO,b),VY(lOO,5),VZ(iOO,b>,T2(100»5) CROS 031 

CROS 03« 
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FIGURE 4. EFFECT OF TRANSITION REGION MODIFICATION ON THE SKIN FRICTION 
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FIGURE 5. TRANSITION CORRELATION CURVE FROM REFERENCE 32 




FIGURE 6. COMPARISON OF EXPERIMENTAL AND CALCULATED TRANSITION LOCATIONS 
FROM VARIOUS METHODS FOR FAVORABLE GRADIENT FLOWS 
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FIGURE 7. FLOW DIAGRAM OF COMPUTER PROGRAM FOR AXISYMMETRIC 
ANALYSIS AND DESIGN METHOD (ADAM) 
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FIGURE 8. SCHEMATIC DIAGRAM OF CYLINDRICAL WAKE SHAPE USED 
TO MODEL SEPARATION 
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FIGURE 9. PRESSURE DISTRIBUTION FOR SPHERE IN SUPERCRITICAL REGION 
USING CYLINDRICAL SEPARATION MODEL 
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FIGURE 12 


COMPARISON OF BASE DRAG CALCULATED BY ADAM TO 
EXPERIMENTAL DATA 



FIGURE 13. SCHEMATIC OF HIGH FINENESS RATIO BODY FROM REFERENCE 35 



FIGURE 14. HIGH FINENESS RATIO BODY AND SIMULATED TUNNEL USED 

IN POTENTIAL FLOW PROGRAM TO ACCOUNT FOR WALL EFFECTS 
ON PRESSURE DISTRIBUTION 
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FIGURE 14. COMPARISON OF CALCULATED "VISCOUS" PRESSURE DISTRIBUTION 
FOR HIGH FINENESS RATIO BODY TO EXPERIMENTAL DATA 
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figure 16. COMPARISON OF CALCULATED "VISCOUS" PRESSURE DISTRIBUTION 
FOR SPHERE IN SUPERCRITICAL REGIME* TO EXPERIMENTAL DATA 
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FIGURE . EFFECT OF "VISCOUS" MODELING ON CALCULATION OF LOCAL SKIN 
FRICTION COEFFICIENT FOR SPHERE IN SUPERCRITICAL REGIME 
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FIGURE 18. COMPARISON OF CALCULATED "VISCOUS" PRESSURE DISTRIBUTION FOR 
SPHERE IN SUBCRITICAL REGIME TO EXPERIMENTAL DATA 
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FIGURE . EFFECT OF "VISCOUS" MODELING ON CALCULATION OF LOCAL SKIN 
FRICTION COEFFICIENT FOR SPHERE IN SUBCRITICAL REGIME 
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